A VARIATIONAL REPRESENTATION OF WEAK SOLUTIONS 
FOR THE PRESSURELESS EULER-POISSON EQUATIONS 



EITAN TADMOR AND DONGMING WEI 

<^^> Abstract. We derive an explicit formula for global weak solutions of the one dimensional 

O^l ■ system of pressure-less Euler-Poisson equations. Our variational formulation is an extension 

(—^ ' of the well-known formula for entropy solutions of the scalar inviscid Burgers' equation: since 

q_) , the characteristics of the Euler-Poisson equations are parabolas, the representation of their 

pL^ ' weak solution takes the form of a "quadratic" version of the celebrated Lax-Oleinik variational 

formula. Three cases are considered, (i) The variational formula recovers the "sticky particle" 
solution in the attractive case; (ii) It represents a repulsive solution which is different than the 
one obtained by the sticky particle construction; and (iii) the result is further extended to the 
-multi- dimensional Euler-Poisson system with radial symmetry. 

Ph 

< 
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1. Introduction 

We study the system of pressureless Euler-Poisson equations 

(1.1a) p t + (pu) x = 0, 

q^'. (Lib) (pu) t + {pu 2 ) x = k P E, E x =p, E{-oo,t) = 0. 

0- ' 

If} ■ Here E is the electric field, k is a given physical constant which signifies the type of the underlying 

l/") . forcing. We distinguish between three different cases depending on the sign on k. 

£\\ [ (i) The attractive case, k < 0. Solutions of (jl.ip always breakdown at a finite time, t = t c , 

O ■ where u x (-,t f t c ) — > — oo. 

(ii) The repulsive case, k > 0. It was shown in [EnLTOl] that in the presence of repulsing 
forcing, there is a large class of so-called sub-critical initial data, u' > —y/2kpo, for which (ll.ljl 
admits global smooth solutions governed by 



(1.2) ut + uu x = kE; 

see [LT021 ILT031 |ChengT08t lTWe08l ICTiaeT08l ILTWIOl WelOal IWelOb] for the critical threshold 



phenomena in related Euler and Euler-Poisson systems. For super-critical data, however, the 
repulsive Euler-Poisson solution breaks-down at a critical time, t = t c < oo, after which (II. 2p is 
no longer equivalent to (ll.lbj) . 

(iii) Finally, there is the neutral case k = 0, governed by the pressureless Euler equations, 
(1.3a) p t + (pu) x = 

(1.3b) (pu) t + {pu 2 ) x = 0. 
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For increasing data, u' > 0, solutions of (jl.3p remains smooth (corresponding to sub-critical 
data in the limiting case k = 0), and are governed by the decoupled inviscid Burgers equation 

(1.4) u t + uu x = 0. 

It is well known that for general non-increasing initial data solutions of (jl.4p will lose their initial 
C 1 regularity at a finite time |La57j . Thereafter, (|1.4p and (|1.3h>[) are not equivalent. Solutions of 
Burgers' equation (|1.4p past the critical time develop shock discontinuities and are given by the 
celebrated Lax-Oleinik formula [Ev98j . Solutions of the pressureless Euler system (|1.3|) develop 
5-shocks. Their construction attracted great attention in the 90's: they were obtained using 
sticky particles formulation in [Z7CH ICPY901 IBG981 ICKR07j . using a generalized variational 
principle in [ERS96J , and through a vanishing pressure limiting process [Bou931 ICL031 ICL04] 
and the references therein. Uniqueness was proved in [HWOlj . 

Here, we are interested in the global weak solution of the (non-neutral) Euler-Poisson system 
(jl.ip . in either the attractive case, k < 0, or in the repulsive case n > subject to general initial 
data, beyond the global regularity in the sub-critical regime studied in [EnLT Ol], Among the 
few known results we mention the existence and uniqueness result of global weak solutions for 
the pressureless Euler-Poisson system [NTu08] . 

In this paper, we construct an explicit formula for weak solutions of the Euler-Poisson (II. ip . 
For the inviscid Burgers' equation, generic shock develops due to the intersection of straight 
characteristics, after which entropy solutions are given by the variational Lax-Oleinik formula. 
For the pressureless Euler-Poisson equations, characteristics become quadratic, and the solution, 
p(x, t), pu(x, t) will be expressed in terms of the minimizer, y{x, t) of the weighted quadratic form 

(1.5) y(x,t) = supjy y = arginf f Q x ,t(s)p (s)ds\, Q x ,t(s) := s + tu (s) + ^-KE (s)t 2 -x. 

y (. y Jo J 1 

Theorem 1.1. [L initial density] Consider the pressureless Euler-Poisson system (jl.ip subject 
to initial data uq(x) := u(x, 0) G C 1 (1R) and < po(x) := p(x,0) G L^R). Set Eq(s) : = 



po(w)dw as the corresponding initial electric field. Then, (p, pu) = (d x R,d x M) is a weak 
solution of (jl.ip . where R(x,t) = R(y(x,t)) and M(x,t) = M(y(x,t)) are given by 

(1.6) R{x,t) = J po{s)ds, M(x,t) = J (u (s) + KE (s)t)p (s)ds. 

Remark 1.1. The presence of Poisson potential is responsible for 'converting' the straight char- 
acteristics familiar from the Burgers' equation, into parabola in Euler-Poisson equations. This 
is reflected in the variational formula M.b]) through the additional term ^KEo(s)t 2 in S1.5\) . 



Remark 1.2. The representation formula (jl.6|) applies for all k£1. When k = 0, one recovers 
the variational formulation of [HWOlj for the pressureless Euler equations. In the particular case 
of k = and po = 1, (|1.5p . ([1.6p become 

( f y (y - x) 2 x 2 1 

y(x, t) := sup <y y = arg inf / tu (s)ds H — L 

y I y Jo 1 1 ) 

recovering the celebrated Lax-Oleinik formula. 

Remark 1.3. For k < 0, our representation formula (|l,6p gives the "sticky particle" solution of 
the pressureless Euler-Poisson system. For k > 0, the weak solution given by l\1.6\) is different 
from the one corresponding to the "sticky particle" model. We discuss the details in sections^ 
and\E 
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We can extend the result of theorem 1 1,1 1 to more general initial measure data, p(x, 0) G ^#+(IR) 
and m(x,0) = (p(x, 0)u(x, 0)) G ^T(R). 

Theorem 1.2. [Measure initial density] Consider the pressureless Euler-F 'oisson system 
subject to initial data, m(x,0) = p(x,0)u(x,0) G ^#(IR) and po(x) :=p(x, 0) G ^+(181), suc/i i/mi 
po(x)dx < oo and tio(^) «s piecewise continuous. The corresponding initial electric field, 
Eq, is given by the average rule 



(1.7) 



E (s) :-- 



1 



p (w)dw + / p (w)dw , 



Then, (p, pu) = (d x R,d x M) is a weak solution of where R = R(y(x,t)) and M = 

M(y(x,t)) are given by 

ry(x,t)+ 

p (s)ds, if Q x ,t{y{xit)) < 0, 



(1.8a) 



and 



(1.8b) 



R{x,t) 



y(x,t)- 



Po(s)ds, if Q x ,t(y(x,t)) > 0, 



M(x,t) 



< JO 



^ JO 



y(x,t)- 



u (s) + KE (s)t)p (s)ds, if Q x ,t(y(x,t)) < 0, 



u (s) + KE (s)t)p Q (s)ds, if Q x ,t(y(x,t)) > 0. 



We note in passing that since y(-,t) is monotonically increasing, consult lemma [2T2l below, the 
one-sided limits in (jl.8j) . y(-,i)±, are well-defined. 

These results can be extended to the weighted multi-dimensional Euler /Euler-Poisson systems 
with symmetry, which is the content of our third main theorem. 

Theorem 1.3. Consider the n-dimensional weighted Euler-Poisson equations 
(1.9a) p t + V • (pu) = 0, u(-, t) : R n H- R n , 

(pu) t + V • (pu u) = Kp|x| n_1 VF, = p, 



(1.9b) 



X 



subject to spherically symmetric initial data, po(x)uo(x) :=po(|x|)uo(|x|)- — -, such that 



x 



< p(x,0) = po(|x|) G L 1 (W l ), and u (\x.\) G C 1 (M + ). Then the system ( TO]) admits a weak 
radial solution, (p, pu), of the form 



(1.10a) 



p(x,i) 



1 dR(r, t) 



. . . , x 1 dM(r,t) 
p(x,t)u(x,t) = -— T - 



<9r 



r := x 



//ere, R(r,t) = R(y(r,t)) and M(r,t) = M{y(r,t)) are given by, 



(1.10b) R(r,t) 



y(r,t) 



po(s)s n ~ 1 ds, M(r,t) 



y{r,t) 



where Eq, is the initial electric field, 
(1.10c) £*,(r):=r n -V r (|x|,0) 



„n-l 



po(s) (w (s) + nE (s)t )s n L ds, 



p (s)ds, r=|x|, 
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and y = y(r, t) is determined in terms of the quadratic form Q r ,t(s) = s + tuo(s) + ^KEo(s)t 2 — r, 

(l.lOd) y(r,t) := sup \y y = arginf f Q r>t (s)s n ~ 1 p (s)ds\. 

y I y Jo J 

Remark 1.4. If n = 0, theorem \1.3\ yields the formula for a global weak solution of the multi- 
dimensional pressureless Euler system with symmetry. If n = 1, it recovers the formula for the 
ID pressureless Euler-Poisson system, that is, theorem \l.l\ with a anti- symmetric Eq. 

The paper is organized as follows. In section [21 we develop the formula for the attractive 
pressureless Euler-Poisson system. In section [3j we extend the formula to the repulsive system, 
and explain different physical meanings of the weak solutions for k < and k > 0. In section [H 
we extend the formula to the multidimensional weighted Euler-Poisson system. 

2. Weak solutions of the attractive pressureless Euler-Poisson system 

2.1. From continuum to particles. We begin with the transport of the center of mass in the 
smooth case. 

Lemma 2.1. Assume that {p{- , t) , pu(- , t)) is a smooth solution of the pressureless Euler-Poisson 
system M-l\) for t £ [0, T]. Let x(a, t) denote the particle path emanating from a£l. Then the 
center of mass of particles emanating from the interval [a, b] is given by the parabola 

(2.1) X(t) = X + U t + ^KEot 2 , t e [0, T] 

where 



rb i-b 

I xpo(x)dx / po(x)uo(x)dx 

J a t j J a 

po(x)dx / po(x)dx 



(2.2) X = ^- b , U = ^— b , E = -( j Po (x)dx+ j Po (x)dx 



1 
2 

Proof. We use the method of characteristics to obtain an explicit solution of (jl.ip . [EnLTOl] . 
Along the particle trajectory, the equations of x and u are 

dxicv., t\ 

(2.3a) — 3^ — = u(x(a,t),t), x(a, 0) = a, 
and 

(2.3b) du{x(a,t),t) = K y^ x ^ t) = kE ( x ^ t ^ t ^ 

where 

(2.4) E(x(a,t),t) = f ^ p(£)de. 

J — oo 

Since 

d_ 
dt 



E(x(a,t),t) = j t (x(a,tj) ■ p(x(a, t), tfj + J ( ^ 



x(a,t) 



u(a, t)p(x{a, t),t)) - / (p(£, t)u{£, t)) d£ = 0, 



the electric field remains constant along x(a,t). Let 

E Q (a) = Po(s)ds, u = u(x(a, 0), 0), p = p(x(a, 0), 0). 
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Then (|2,3b|) can be simplified to 

du(x(a,t),t) 

= KEc\, 

dt ' 

which yields 

(2.5) u(x(a, t),t) = uq + nE^t. 
This together with (|2.3ap yield the equation of the particle path 

kEq^ 

(2.6) x(a,t) = a + uot -\ — . 

Let 

tv dx , K Po t 2 , du {a) 

Taking the x derivative of u(x(a,t),t) yields 

(2.7) „, Wa , ,), ,) = = <±™< _ 

oa oa , Kpnt 1 (a, t) 

l + u' t+-^- 

Plugging (|2.7p into the mass equation (jl.lap yields 

d , TAa,t) 
—p(x(a,t),t) = -u x p = -— — —p(x[a,t),t). 
dt 1 (a, t) 

Solving this equation, we obtain 

Po 



(2.8) p(x(a,t),t) - 

T(a,t) 

Equipped with (12. 5p . (|2.6I) and (j2.8|) . we can find the position X(t) of the center of gravity of 
the mass on [x(a,t),x(b,t)\ at time t 

rx(b,t) „{, 

£p(£,t)d£ / x(a,t)p(x(a,t),t)T(a,t)da 



X(t) 



x(a,t) 



X(b,t) ft^t E (b)-E (a) 

x(a,t) 

f b ( t \ KE (a)t 2 \ . , , „ f b „ . . . . , 

/ ^a + uo(a)tH J/Oo(Q:)aa «i / £/o(a)po (ajda 

~ E (b)-E (a) =X + U t+ 2{ E o{b) _ E(j{a)) 



,6 

Kt 2 / ^ (a)d-Eo(a) 

J a 



Xo + C/ ° t+ 2(£b(6)-Eb(a)) - Xo + Uot+ 4(E (b) - E (a)) 
X + U t + ^KE t 2 . 



□ 

Lemma [2. II tells us that if we replace the initial mass along [a, b] by a Dirac mass, J po(x)dx, 

situated at location Xo, give it velocity Uo, and apply to this particle an electrical field Eq, then 
the trajectory of this particle is the same as the trajectory of the center of mass of the whole 
initial interval [a, 6]. Motivated by this lemma, we continue to deal with measure densities. In 
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particular, if the density p(-,t) at location x has a Dirac mass with strength m(x), then we set 
the electric field at that point to be 



E(x,t) := I p(s,t)ds + 



m(x) 



/X 
p(s,t)ds. Combining both cases, we arrive at the following definition 
-oo 

of the electric field, which is in agreement with (|1.7p . 

Definition 2.1. [The electric field] We define the electric field, E(x,t), as 

(2.9) E(x,t):=-(j p(s,t)ds + J p(s,t)dsj, 

Equipped with definition 12. 1\ we will show later that the conclusion of Lemma 12.11 remains 
valid even after collision takes place. Therefore, setting E as in (12.9|) guarantees that the 
trajectory of the center of mass is independent of whether there is a collision or not. This is a 
key point which explains the validity behind our approach. 

2.2. The dynamics of two Dirac masses. Given the density p(-,t) as a non-negative measure 
and the velocity u(-,t) which is uniformly bounded, we introduce the corresponding mass and 
momentum of the system, which play an important role throughout the paper, 

/X + I-X + 

p{s,t)ds, M(x,t) := / p(s,i)u(s,t)ds 
-oo J — oo 

To illustrate the construction of a weak solution solely from the physical principles, we start 
with the simplest example of two particles governed by an attractive force, k < 0. We con- 
sider two particles with masses mx,m2 at initial positions y±(0) < 2/2(0), and initial velocity 
iti(O), 1*2(0), respectively. Thus, the initial density and momentum consist of two Dirac masses 

po(x) = mi5(x - yi(0)) + m 2 S(x - 2/2(0)), 
p (x)u (x) = miui(0)5(x - yi(0)) + m 2 U2(0)5(x - y 2 (0)). 

If there is no collision, the electrical fields which are associated with the first and second particles 
are, respectively, E\ := ^ and E 2 := m\ + Hence the velocity and location of the first 
particle are 

(2.11) «i(*) = u 1 (0) + K^t, y l (t) = y 1 (0)+u l (0)t + ^ 1 ^-t 2 . 
The velocity and location of the second particle are 

(2.12) u 2 (t) = u 2 (0) + k (mi + ^) t, y 2 {t) = y 2 {0) + u 2 (0)t + -k (mi + ^) t 2 , 

Therefore the velocity and location of the center of gravity of the system are 

, . miu\{t) + m 2 u 2 (t) _ miUi(O) + m 2 u 2 (0) 1 m\ + 2mim 2 + m\ 

U[t J — — "T" — K 1 

m\ + m 2 mi + m 2 2 m\ + m 2 

m\ , mi+m 2 miui(0) + m 2 -u 2 (0) 

= u(0)+k 1, u(0) := ■ , 

2 mi + m 2 



and 

y{t) 



WEAK SOLUTIONS OF PRESSURELESS EULER-POISSON EQUATIONS 



miyi(t) + m 2 y 2 (t) 



mi + 7712 

toi2/i(0) + m23/2(0) , mitti(O) + m 2 u 2 (0) 1 mf + 2mim 2 + m| 2 

mi + m 2 mi + m 2 4 mi + m 2 

/ n x . / n x. . 1 "ii + m 2 ,2 / n x mij/i(0) +m 2 j/ 2 (0) 
= y(0)+n(0)t+- K ^^t , y(0):= . 

If collision takes place at time t, then a new particle will be formed at y(t) with velocity u(t): 
the motion of this new particle is determined by 

(2.13a) u new (t) = u(t) + k (t-t) = u(t), 

and 

(2.13b) y new (t) = y(t) + u{t)(t + l K mi + m V = y(t). 

Hence, when and where the collision occurs, it will not change the trajectory of the center of 
mass of this system. 

The standard way to describe the state of this two-particle system at any time t, would be to 
provide the following information: 

(i) if the two particles have not collided before time t: provide the position and velocity of each 
particle at time t; 

(ii) if the two particles collided at some time before t: then they are "stuck" as one particle 
thereafter, and we provide the position and velocity of this new particle at time t. 

Alternatively, we can provide a complete description of the two-particle system in terms of 
the mass R(x,t) and momentum M(x,t). To this end, take the difference of (|2.12p and (|2.1ip . 
to find 

y 2 (t) - yi(t) = \k + ^) t 2 + mo) - «i(0))< + (y 2 (o) - jfc(o)). 

Since the two particle system is attractive, k < 0, there is a positive t c > 0, such that, 

y 2 {t) > yi(t) 0<t<t c , 

V2(t) = yi(t), t = t c , 

y 2 (t) <yi(t), t> t c . 

Then, there are four possibilities to consider, depending on the relative positions of qj(t) := 
yj + ujt + ^KEjt 2 , j = 1, 2: 

(i) if both qi(t) and q 2 (t) are to the left of x : qi(t),q 2 (t) < x, then R(x,t) = mi + m 2 and 
M(x, t) = m\{u\ + nEit) + m 2 {u 2 + nE 2 t); 

(ii) if both qi(t) and q 2 (t) are to the right x : qi(t),q 2 (t) > x, then R(x, t) = 0, and M(x, t) = 0; 

(iii) if qi(t) < x and q 2 (t) > x, then R(x,t) = mi, and M{x,t) = mi{ui + nEit); 

(iv) finally, if qi(t) > x and q 2 {t) < x, this means that collision occurred earlier, and there is a 
new particle with mass m = mi + m 2 . The velocity and location of this new particle at time t 
are given by (|2.13a|) . (|2.13b|) . 

miiux + nExt) + m 2 (u 2 + nE 2 t) miqAt) + m 2 q 2 (t) 

«W = : . V{t) = • 

mi + m 2 mi + m 2 
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Therefore, in this case R(x,t) and M(x,t) are determined by the sign of the expression 
TOiqi(t) + m 2 q 2 (t) _ _ m 1 (q 1 (t) - x) + m 2 (q 2 (t) ~ g) 
mi + m 2 mi + m 2 

namely 

J?(z, t) = mi + m 2 | if ~ + m 2 (g 2 (*) - x) < 

M(x, t) = mi(«i + KE x t) + m 2 (u 2 + nE 2 t) J ' m i + m 2 

R(x,t) = M(x,t) = 0, otherwise. 
Summarizing the four cases above we observe: 

(i) if qi(t) < x, then no matter whether collision happened or not, mi will contribute to the 
mass R(x, t); 

(ii) if qx(t) > x then the contribution of mi to the mass R(x,t) depends on whether the second 
particle is slow enough, that is, whether mi(qi(t) — x) + m 2 {q 2 {t) — x) < 0. 

2.3. The dynamics of general mass distribution. In this subsection, we extend the defini- 
tion of R and M to the general attractive case. Suppose that there exists Y\ such that 

q 8 (t) := s + u (s)t + ^KE(s,0)t 2 < x, Vs < Y x . 
Then, independently whether collision occurred or not, the particles emanating from y < Y\ will 
end up to the left side of x at time t, and therefore, the part of the mass / po(s)ds will be 

J ~ oo 

on the left side of x at time t. If, on the other hand, 

q s (t) = s + u (s)t + ^KE(s,0)t 2 > x, Vs G (Yi,y 2 ], 

then the position of the mass of this part, / po(s)ds, relative to x, depends on whether there 

is enough slow material which will collide with this part. That is, if there exists Yjj > Y 2 , such 
that 

,y 2 nY 3 nY 3 

s (t)p (s)ds+ / q s (t)p (s)ds / q s (t)p (s)ds 
(2.14) = hi < Xi 

/ p (s)ds+ / po(s)ds / po{s)ds 

JYx Jy 2 Jy 1 

then the mass on (Yi, Y 2 ] will be end up to be on the left side of x. We can rewrite (|2.14p in the 
equivalent form 

Qx,t( s )Po(s)ds < 0, Q x ,t{s) = q s (t) - x. 

rY 2 

This is similar to the case (iv) of the two particle system, where / po(s)ds corresponds to the 

Jy 



r*3 

first particle, and / po(s)ds corresponds to the second particle. Continuing this process, we 
JY 2 

find that the exact amount of mass which ends on the left side of x, is given by 

rvfat) 

(2.15) R(x,t) := / p (s)ds, 



oo 
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where y(x,t) is determined by (jl.5p . 

(2.16) y{x,t) = sup ly y = arginf f Q x ,t(s)p (s)ds\, Q x ,t( s ) = Qs(t) - x. 

v I y Jo J 

The momentum then follows 

po(s)(«o(s) + W 

-oo ^ ' 

2.4. Proof of theorem II. II for the attractive case. As a preparation for the proof of theorem 
ll.ll with k < 0, we first characterize the entropy solution of (jl.ip in terms of a one-sided Lipschiz 
condition. This is the content of the following two lemmas. 

Lemma 2.2. Let y(x,t) be the minimizer in il.5\) . Then y(-,t) is non- decreasing 

y(xi,t) < y(x2,t), Vxi < xi. 
In 'particular, y(x+,t) := lim z ^ x+ y(z, t) is well-defined and y(x,t) = y(x+,t), Vx. 
Proof. Assume there exists x\ < X2 such that y(x\,t) > y(x2,t). Since y(xi,t) minimizes 



/ Qx!,t{s)po{s)ds, we have / Qx 1 ,t(s)po(s)ds < 0. It follows that 

JO Jy(x 2 ,t) 

ry(xi,t) ry(xi,t) ry(^i,t) 

/ Qx 2 ,t( s )Po( s ) ds = / Qxi,ti s )Po( s )d s + / ( x i ~ x 2 )po(s)ds < 0. 

Jy(x 2 ,t) Jy(x 2 ,t) Jy{x 2 ,t) 

That is, 

rvipit) ry(a>2,t) 

Qx 2 ,t(s)po(s)ds < / Q X2t t(s)p (s)ds, y{x±,t) > y(x 2 ,t), 



I 

Jo 
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which is a contradiction to the definition of y(x2,t). 

ry(x+,t) 

Next, assume there exists x such that y(x+,t) > y(x,t). Then / Qx,t(s)po{s)ds > 0. 

Jy(x,t) 

Therefore, for e > small enough, we have 

ry(x+e,t) 

Qx+e,t( S )P0(s)ds 

y(x,t) 

y(x+,t) ry(x+e,t) 



ry(x+,t) ryyx+e,t) 

I Q x+€)t (s)po(s)ds+ / Q x+et t(s)p (s)ds 

Jy(x,t) Jy(x+,t) 
ry(x+,t) ry(x+,t) ry(x+e,t) 

/ Qx,t(s)po(s)ds - / e(s)p (s)ds+ / Q x+e j(s)p (s)ds > 0. 

Jv(x.t) Jv(x.t) Jvix+.t) 



That is 



y(x,t) ry(x+e,t) 

Qx+e,t( S )P0(s)ds < / Q x +e,t( S )P0(s)ds, 



which is a contradiction to the definition of y{x + e, t). □ 

Lemma 2.3. [One sided Lipschiz condition] Consider the attractive Euler-P oisson system (jl.ip 
with K < 0. Ify(x,to) > y(x—,to), then the values in the open interval 7 G (y(x—,to),y(x,to)), 
cannot be reached by evolving (jl.ip along particle path, namely, 7 / a minimizer y(x,t). 

Thus, if y(x,to) > y(x—,to), then according to (jl.6p . p(x,to) will be a Dirac mass and the 
lemma [2731 tells us that once a Dirac mass is formed, it will never split. 
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Proof. By the definition ()1 . 5[) . to prove the Lemma, it is enough to show V7 G (y(x— , to), y(x, to) 

pw /"J 

Vz G R, Vt > to, there exists to such that / Q z ,t(s)ds < / Q z ,t(s)ds, that is, 7 does not 

jo Jo 

minimize / Q z ,t(s)po(s)ds. For every 7 G ( j/(a;— , to)? y(x, to) ) > we have 







(2.18) / Qx,t (s)Po(s)d«<0< / Qx,t (s)Po(s)ds, 

•Ay Jy(x-,t ) 

These inequalities can be rewritten in terms of the quadratics /i(t) and /2(t), 

/ 2 (*o) < a < /i(t ), fj(t) := a 3 t 2 + bjt + c,-, 

/"Y 

where the coefficients of fj(t)'s are given by in terms of m\ := / po(s)ds and m 2 

J ,in) 



y(a;,to) 

p (s)ds: 



1 



1 P 1 IP IP 

ax = — / -nE (s)po(s)ds, h = — / u (s)p (s)ds, c x = — / sp (s)ds, 



m l J y(ai- to) 2 m l Jy(x-,t ) mi Jy( x - to ) 

and 

I ry{x,t ) 1 1 ry{x,t ) \ rv(x,to) 

a 2 = — / -KE (s)p (s)ds, b 2 = — / u (s)p (s)ds, c 2 = — / sp (s)ds. 



m 2 J 1 2 m 2 J 7 m 2 j~ f 

Notice that a 2 < kEq{^) < a\ < 0, and c 2 > 7 > c\. Setting /i(t) = fi(£), we find two solutions, 
one positive and one negative. Denote the positive one by t c , then 

/l(t)</2(*)> 0<t<t c 

/i(t) > /2W, * > * c . 

It follows that to > t c . For every t > to and z G R we have 

7 

Q z ,t( s )Po(s)ds = mi(ait 2 + M + ci - z) = mi(/i(t) - z), 

2/(a-,to) 
/•v{x,to) 

/ Qz,t(s)p (s)ds = m 2 (a 2 t 2 + 6 2 t + c 2 - z) = m 2 (f 2 (t) - z). 

If /i(t) > z, then 

ry(x-,t ) 



I Qz,t{s)po{s)ds > / Qz,t( s )Po(s)ds; 
Jo Jo 



if fl(t) — then f 2 (t) < /i(t) < z, which implies 

Qz,t{s)po{s)ds > / Q Zji (s) / o (s)c?s- 
'0 Jo 
ry 

Therefore 7 does not minimize / Q Z! t{s)po(s)ds. □ 

Jo 

Remark 2.1. It zs straightforward to verify that for smooth initial data, uq(x) G C(R), we have 
Qx,t(y(x,t)) = 0, i.e., 

y(x, t) + tu (y(x, t)) + ^E (y(x, t))t 2 = x, Vx. 
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Remark 2.2. Consider two adjacent discontinuous points a and b of y(x,t), i.e., y(a—,t) < 
y(a,t), y(b—,t) < y(b,t) and y(x—,t) = y(x,t) for every x G (a, b). Then R and M are 
continuous on (a, b) at time t. Moreover, combining lemma W^ and remark [2J[ we find that the 

characteristics emanate from the interval I y(a, t),y(b, t)j at t = 0, 



x(t) = x + u (x )t + ^KE (x )t 2 , Vx G (y(a, t),y(b, t)), 
will not intersect before t. 

As a final preparation for the proof of theorem 11.11 we calculate the distributional derivative 
of jump discontinuities across curves over surfaces, which will be useful when dealing the singular 
part of Rt and M x . We summarize this calculation in the following lemma. 

Lemma 2.4. Consider an open region, V C R 2 , and a curve, C(t,x(t)) in the t — x plane which 
divides V into two parts, V 1 and V r , and assume that a function S(t,x) is smooth on either 
side of this curve C, with values S l on V 1 and S r on V r . The weak derivative of S is given by, 

S x = S l x + S r x + (S r - S l )v2VC s t = s l t + S r t + (S r - S l )u lVc . 

Here, v is the outward unit normal vector of V 1 on boundary C, 

= ( ) - ( "I ^ 

{Ul,V2h ^ 2 (*) + iV(* 2 (*) + V " 

and rjc(t,x) is a surface measure supported on the curve C, satisfying 



<f>(t,x)ric(t,x)dtdx = / <t>(t,x(t))y/x 2 (t) + Idt, \/(f> G C^°(IR 2 ). 

R 2 J a 

Proof. For every <j)(x,t) £ C^°(y), 

4> x (x,t)S(x,t)dtdx = / (p x (x,t)S l (x,t)dtdx + / (p x (x,t)S r (x,t)dtdx 
V Jv 1 Jv r 

x,t)S l x (x,t)dtdx + / 4>(t,x(t))S l (t,x(t))u 2 dC - / <j>{x,t)S r x {x,t)dtdx 
v l Jc Jv r 



<f>(t, x(t))S r (t, x{t))v 2 dC (where (dC = ^(dx/dt) 2 + 1)) 

0(x,t)S l x (x,t)dtdx - [ <j)(x,t)S r x (x,t)dtdx - [ 0(t,x(t))(s r (t,x(t)) - S l (t,x(t)))u 2 dC, 



JV r JC 



<j>(x,t)S l x (x,t)dtdx - / (p(x,t)S r x (x,t)dtdx - / (f>(t, x(t)) ( S r (t, x(t)) - S l (t, x(t)) ) v 2 r) C dtdx. 
v l Jv r Jv ^ 

Hence S x = S l x + S x + (S r — S l )v 2 rjc- Similarly, 

(f>t(x,t)S(x,t)dtdx = / 4>t(x,t)S l (x,t)dtdx + l <f) t (x,t)S r (x,t)dtdx 



(j)(x,t)S l t (x,t)dtdx + / (t,x(i))MC - / (f>(x,t)S r t (x,t)dtdx 

v l Jc Jv r 

0(t,x(t))5 r (t,2;(t))i/ 1 dC 

<p(x,t)S l t (x,t)dtdx - / <p(x,t)Sl(x,t)dtdx - / <t>(t,x(t))(s r (t,x(t))-S l (t,x(t)))v 1 dC, 
v l Jv r Jc ^ ' 

<f)(x,t)S l t (x,t)dtdx - / 4>(x,t)Sl(x,t)dtdx - / (j)(t,x(t))(s r (t,x(t)) - S l (t,x(t)))uii] C dtdx. 
v l Jv r Jv ' 
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Hence S t = S l t + 5 t r + (S r - S l )vir) C . □ 

Equipped with lemmas 12.31 and 12. 4| we are now ready to complete the proof of theorem 11.11 
in the attractive case, k < 0. 

Proof of theorem li.il with k < 0. Step #1 [the mass equation]. First, we show p and pu satisfy 
the mass equation (jl.lap in the weak sense. We need to verify that 

poo poo poo poo poo poo 

(2.19) / / pt4>+(pu) x (f)dxdt = / / R xt (f> + M xx cf)dxdt = / / R(f> xt + Mcf) xx dxdt 

JO J-oo J J-oo JO J-oo 

vanishes for all test functions (j) £ C^°(\0, oo) x (— oo, oo)^ . To this end, we claim that 

dR dM 

< 2 ' 20 > «— 5T 

It follows (|2.20p that there exists a function G such that 

\ dG(x,t) , ,. . dG(x,t) 

(2.21) fl(ac,t) = M(a,t) = K m ! . 

By plugging (12.211) into (|2.19p we obtain 

/"OO f'QO f'OO POO , v 

(2.22) / / Pt <t> + {pu) x( t)dxdt= / / G0 xxt - G<f) xxt dxdt = 0, G C c °°(lR x K+). 

We decompose -Rj and M x into the sum of an absolutely continuous measure and a singular 
measure. We denote 

R t =:Rf + Rf, M X =:M« + M S X , 

where Rf and M x are absolutely continuous measures, Rf and M x are singular measures. More 
precisely, if (x,t) is a jumping point (that is, y(x—,t) < y(x,t)^, then Rt and M x are singular 
at (x,t); if y(x—,t) = y(x,t), then i? and M are continuous, hence i?j and M x are not singular. 
The detailed proof of (I2.20p is carried out by verifying that Rf = —M x and R s t = —M x in steps 
#l(a) and #l(b) below. 

Step #l(a) [the regular part of the mass equation]. We show that Rf = —M x . If (x,t) is 
a jumping point, then Rf = —M x = 0. Otherwise, y(x,t) is continuous at (x,t), and by the 
definition of R and M, we have 

and M? = on(v) (un(v) + KEn(v)t) ^ 



(2.23) Rf = Po(y)~^, and M a x = p (y)(u Q (y) + KE (y)t) g 

By remark [2T2l the equality Q x ,t(y(x, t)) = is valid in the neighborhood of (x, t) : Q x ,t(y(%i t)) 
tuo{y) + y + -nEo(y)t 2 — x = 0. Taking partial derivatives with respect to t and x yields 

Hence 

9y u {y) + nE (y)t dy 1 



5i d y Q x j{y) dx d y Q x j{yY 

and therefore 

(2-24) ^ = -(n (y) + Kj B (y)t)^. 

Combining (f!T24l) and (f!T23"l) we obtain Rf(x,t) = -Mg(x,t). 
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Step #l(b) [the singular part of the mass equation]. We show that i?| = — MJ. If y(z—,t) = 
y(z,t), then Rf = —M£ = at (z,t). Otherwise, (z,t) is a jumping point where we have 

ry(z,t) / l \ 

(2.25) / (s + tu (s) + -KE {s)t 2 -z )p o (s)ds = 0, 

Jy(z-,t) V 2 J 

(2.26) y(z, t) + tu (y(z, t)) + ^E (y(z, t))t 2 - z = 0, 

(2.27) y(z-,t) + tu (y(z-, t)) + - K S (!/(«-, t))t 2 - z = 0. 



Denote the trajectory of this jumping point by C(t, z(t)). According to Lemma 12.41 we have 

(2.28) R!(z,t) = B(z,t)r ]c ^=M= 
and 

(2.29) M s x {z, t) = -B(z, t)u(z, t)r, c 1 ; 

V 1 + z (*) 



where 



tf (, jt ) / Po(s)(^n (s) + KE (s)tjds 



y(z-,t) 



B(z,t) = / po(s)ds, u(z,t) 

Jy(z-,t) B{z,t) 

Thus, to prove Rf = — at the jumping point (z,t), it remains to show that the propagation 
speed of this jumping point, z(t) =: u(z,t) is actually given by u(z,t). We provide the details 
of u(z,t) = u(z,t) below. 

The location of this jumping point at time t + At is (t + At, z(t + At). Hence 

py(z(t+At),t+At) 1 

(2.30) / (s + {t + At)u (s) + -KE (s)(t + At) 2 -z(t + At) )p (s)ds = 0. 

Jy(z(t+At)-,t+At) V 2 / 

Combining, 

y(z{t + At), t + At) = y(z, t) + O(At), 

(2.31) y(z(t + At)-, t + At) = y(z-,t) + O(At), 

z(t + At) = z + u(z, t)At + 0(At 2 ). 
with (f!T30D and (12351) . we obtain 

r-y(z(t+At),t+At) ^ 



= / s + (t + At)u (s) + -K£ (s)(i + At) 2 -z(t + At) po (s)ds 

Jy(z(t+At)-,t+At) V 2 / 

= / (s + (t + At)u (s) + -KE (s)(t + At) 2 - z-u{z,t)At + 0{At 2 )) )p (s)ds 

Jy(z-,t) V 2 / 

+ / [s + {t + At)uo(s) + -KE (s)(t + At) 2 -z-u(z,t)At))p (s)ds 

Jy(z(t+At)-,t+At) V 2 / 

py{z(t+At),t+At) 1 

+ / s + (t + At)uo(s) + -K J E (s)(t + At) 2 -z-n(z,t)At) /) (s)ds. 
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For every s G ( y(z(t + At)—, t + At),y(z—, t) I , using Taylor expansion and (|2.26j) . we obtain 



8 + tuo(s) + ^E (s)t 2 -z = y{z-,t) + tu (y{z-,t)) + ^ K E {y(z-,t))t 2 -z + O(At) = O(At). 



Hence 

rv(z-,t) , i s 

/ [s + (t + At)ito(s) + -nE (s)(t + At) 2 - z - u(z, t)At) p {s)ds = o(At). 

Jy(z(t+At)-,t+At) ^ 



Similarly 

»y(z(t+At),t+At) 



/ ( s + (t + At)uo(s) + -KE (s)(t + At) 2 - z- u(z, t)At) ) p Q (s)ds = o(At). 



Therefore 

ry(z(t+At),t+At) 



= / s + (t + At)u (s) + -KE (s)(t + At) 2 -z(t + At) p (s)ds 

Jy(z(t+At)-,t+At) V 2 / 

= / s + (* + At)uo(s) + -KE (s)(t + At) 2 - z - u(z,t)At) )p (s)ds + o(At) 

/ 1 \ 

= / (s + tu (s) + -nE {s)t 2 - z))p (s)ds ( — which is by $T25§) 



+ / u (s)At + -K^ (s)(2tAt + At 2 ) - u(z, t)At) p (s)ds + o(At) 

Jy(z-,t) V 2 / 

At / f«o(s) + nE (s)t - u(z,tj)p (s)ds + o(At). 



This concludes our argument that u(z,t) = u(z,t), 

ry(z,t) 



Po(s)[u (s) + kE (s)t)ds 
p (s)ds 



u(z,t) = ^ = u(z,t), 



2/0-*) 

and Rf = — M x follows. Thus, (I2.20P holds, and the mass equation (jl.lap is satisfied in the weak 
sense (12^221 . 

Step #2 [the momentum equation]. Next, we verify the momentum equation (|l.lb|) in a 
similar way. We apply test functions cf> £ C^°( [0,oo) x (—00,00)) to (|l.lb|) . then 



noo 
{pu) t <t) + (pu 2 ) x (p - npE(j)dxdt 
-00 

noo poo roo 

M xt( j) + W xx <p - Z xx cj)dxdt = / / M(j) xt + (W - Z)cf> xx dxdt, 
-00 Jo J-oo 



where 

"X n/r2 



f x M 2 (s t) 

W{x,t):= p 7 \( ds, Z xx (x,t) = k P E. 
Jo R s {s,t) 
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We will show below, in steps #2(a)-2(b), that 

dM d(W - Z) 



(Z33 > m 

This yields the existence of such that 



dx 



(2.34) 



M 



W -Z 



d^ 

dx dt ' 

which in turn, implies that the momentum equation, (jl.lbp . holds in its weak formulation (|2,32p . 

Step #2 (a). The main claim here is that Z x is given by Jq Kpo(s)Eo(s)ds, that is, if we let 
F denote Z r then 



(2.35) 



F x (x,t) = Kp(x,t)E(x,t), F(x,t) :- 



y(x,t) 



Kp (s)E (s)ds. 



As before, we distinguish between two cases. In the case x is a continuity point, y(x—,t) 
y(x,i), then 



(2.36a) 



dF dy 
a- = Kpo{y)E {y)—. 
ox ox 



Combining p(x,t) = —R(x,t) and (11.61) . we obtain 

ox 



(2.36b) 
and 

(2.36c) 



dx dx 



p (s)ds 



dR(x, t) 
dx 



oo 

y(x,t) 



p (s)ds 



oo 
oo 



Po(s)ds 



= p(x,t), 
po(s)ds^J dw 

y(x,t) 



p (s)ds = E (y). 



Plugging (I2.36bp and (|2.36cl) into (I2.36a|) . we obtain ([235]) . F x = Kp(x,t)E(x,t). 

Next, we consider the case of a jump discontinuity at x = z, where y(z,t) > y(z—,t). Then 



dF d /■»(*'*> 

dx x=z dx ./.oo 

= 8{x — z) 

On the other hand 

kp(z,t)E(z,t) = 



Kp (s)E (s)ds 



S(x — z) 



!/(*,*) 



np (s)E (s)ds 



y(z-,t) 



y(z-,t) 



K E (s)dE (s) = 5{x - zY-( E (y(z,t)) 2 - E (y(z-,t)) 2 ). 



dx 
-S(x - z) 



5{x — z) 



2/0-*) 
2/0-*) 



p {s)ds 



Po(s)ds 



p(s,t)ds+ / p(s,t)ds 



dR(s,t) 



-oo 



ds 



ds + 



dR(s,t) 



ds 



ds 



Po(s)ds + 



8(x - zY- ( E (y(z, t)) 2 - E (y(z-,t)) 2 ) . 



y(z,t) 



Po(s)ds 



S(x - z) E (y(z, t)) - E (y(z-,t)) E (y(z-,t)) + E {y(z, t)) 
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Thus, recover (|2.35p F x = Kp(x,t)E(x,t) also at any jumping point z, 

(2.37) Z X = F= /' K p (s)E (s)ds. 

Jo 

Step #2(b). To prove (|2.33p . we decompose M t and (W — Z) x into the sum of an absolutely 
continuous measure and a singular measure, denoting 

M t = M t a + M t s , (W - Z) x = {W- Z) a x + (W- Z)%, 

where {} a denotes an absolutely continuous measure and {} s denotes a singular measure. The 
absolutely continuous measure of (W — Z) x is given by (W x — F), where 

(2.38) ^ := iT ) '°o(s)(wo(s) + ^o(s)t) 2 ^- 

Since y(x, t) is continuous almost everywhere, we have W x = W x almost everywhere. 

Consider the absolutely continuous parts: if y(x—,t) < y(x,t), then M^(x,t) = —{W — 
Z) x (x,t) = 0; otherwise, if y(x—,t) = y(x,t), then 



(2.39a) 



/ \dy r y( - x '^ 

M t a {x,t) = po(y)[u (y) + KE {y)t)-^+ Kp (s)E (s)ds 



Po(y)(u (y) + KE (y)t)^ + F, 



and 

(2.39b) (W - Z) a x (x,t) = W a x -F = p (y)(u (y) + KE (y)tf^ _ p 

Combining (|2.39a|) . (|2.39bl) and (l23i]) . we obtain Mf = -(W - Z)%. 

For the singular parts: if y(x—,t) = y(x,t), then M f s (x,t) = — (W — Z) x (x,t) = 0; otherwise, at 
any jumping point (z,t) we have, (by the definition of R, M and Z), 

M t s (z,t) = B(z,t)u(z,t)7 lc * ( *?, ._ , 

V 1 + z ( t ) 

(W - Z) x (z,t) = W°(z,t) = -B(z,t)u 2 (x,t) Vc - 



When deriving (W—Z) x , we used the fact that (i) Z x = (since Z x = F is bounded everywhere), 
and (ii) we have already shown, z = u(z,t). Hence, we also have for the singular part M t s = 
— (W — Z) s x . Thus, the momentum equation (|l.lbp holds in the weak sense. □ 

The proof of theorem 11.11 covers of course the case of smooth solutions and in particular, the 
variational formula (jl.6j) describes the globally-in-time smooth solutions in the sub-critical case 
[EnLTOl] . In this sense, (|1.6h could be viewed as an extension that covers both sub-critical and 
super-critical initial configurations. We close this section by reproducing the proof of theorem 
l.ll for the simpler case of smooth solutions. 



Lemma 2.5. Consider with a smooth solution on [0, T). Then (p, pu) given in U.6\) is 

that smooth solution. 

Proof. We denote the smooth solution by (p,u), the solution given by (jl.6p by (p,u), and we 
show (p(-,t),u(-,t)) = (p(- , t) , -u(- , t)) , t £ [0, T). We can solve the equation by the method of 
characteristics, and we have shown the details in the proof of Lemma 12. 11 Since the solution is 
smooth, no characteristics will intersect. For every given x G R, t G [0,T), there is an unique 
y(x,t), such that the characteristic which emanates from y(x,t) arrives x at time t. We claim 
y(x,t) is the unique minimizer of Jq Q x ,t(s)po(s)ds, that is, 
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y(x,t) = arginf / Q Xtt (s)p (s)ds. 
y Jo 



Otherwise there exists w ^ y such that 

ry(x,t) 



Qx,t( s )Po( s ) ds - / Qx,t( s )Po(s)ds 



i.e., 



/ Qx,t{s)po{s)ds > 0, if w < y(x,i), 

J w 

Qx,t( s )Po(s)ds < 0, if w > y(x,t). 



Iy(x,t) 

Combining the above inequalities with Q x ,t(y) = 0, we obtain: if w < y, then 3z £ (w,y), 
such that Q x ,t(z) > 0, which implies the characteristics emanate form z and y must intersect 
no latter than t; otherwise w > y, then 3z 6 (y,w), such that Q x ,t{ z ) < 0, which implies the 
characteristics emanate form z and y must intersect no later than t. This is a contradiction to 
smooth solution exists up to T > t. Therefore, 

rvixt) ry(x,t) 

R(x,t) := J Po{s)ds, M(x,t) := J \u (s) + KE (s)tj p (s)ds, 

and 

dy(x,t*) ^ dy(x^t^) ~ ( \ 

p{x, t) = — g^—po(y)> p( x > f ) = — d^~ p °( yS ) { u °( yS ) + ^Eoiyjt) . 

Thus 

u(x,t) = u (y) + nE (y)t = u(x,t). 
To verify p(x,t) = p(x,t), it is enough show that 

p(s,t)ds= / Po{s)ds( = R(x,t) = / p{s,t)ds), Vx, Vt < T. 

-oo >/ — oo V J —oo J 

We denote the inverse function of y by x(y,t), which means the characteristic starts from y 
arrives x at time t. Then (|2.40p is equivalent to 

rx(y,t) 



rx{y,t) rv 
(2.41) / p(s,t)ds= / p (s)ds, My, \ft<T. 

J —oo J —oo 



Notice that y = x(y,0), so (|2,4ip equivalents to: for every fixed y, E(y,t) := /^^' p(s,t)ds is a 
constant. Physically, this is clear: since no charge can across the particle path, by conservation 
of charge, E remains constant along x(y,t). To show E is a constant along x(y,t), we take the 
time derivative of E, then 



0/ 



rz(?/>*) 

t)p(x, t) — I ( p(s, i^x, i) ) ds = up — pu = 0. 



□ 
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3. Weak solutions of the repulsive pressureless Euler-Poisson system 

For k > 0, (jl.6p still yields a weak solution of the Euler-Poisson system, although it may 
be different from the one corresponding to the "sticky particle model". The following example 
demonstrates this point. 

Example 3.1. Let k = 1. We consider a system of two initial Dirac masses: both have mass 
1, their initial positions are and 1, their initial velocities are 2 and 0, respectively. Hence 
the initial density and momentum are po(x) = 5(x) + 5{x — 1), and po(x)uo(x) = 25{x). The 
characteristics of them are 

y l (t) = 2t+h 2 , y 2 ( t ) = l + h 2 . 

If we set yi(t) = y 2 {t), then it yields two solutions: t\ = 2 — y/2 and t2 = 2 + \/2. The dynamic 
of the weak solution given by our formula is: the two particles collides at time t\, then they stick 
as one particle, when this new particle arrives the location yi(t 2 ) = 2/2(^2) at time t 2 , it splits 
into two particles again. 

Therefore, in the repulsive case k > 0, the weak solution given by (|1.6jl subject to smooth 
initial data — even if the data is super-critical, will eventually become smooth again; consult 
remark 13.21 below. 

To further clarify the different behavior of solutions to the attractive and repulsive pressureless 
Euler-Poisson, consider a two particles Euler-Poisson system with mass mi,rri2, initial velocity 
U\,u 2 , initial position y\ < y 2 . If K < 0, then 

1 mi 9 1 !Ti2 \ 9 

(3.1) yi{t) = Vi + Uit + -K—t = y 2 (t) = y 2 + u 2 t + -«(mi + — )t 

always yields one and only one positive solution t c , such that 

l/i (t) < lfe(t), < t < t e 

yi(t) > jfe(t), t > t c . 

If k > 0, however, then (|3.ip can either have no positive solution or two positive solutions, and 
for this is the reason, our representation formula gives a "non-sticky particle" solution for the 
repulsive model. 

We now turn to the proof of theorem 11.11 in the repulsive case, k > 0. 

Proof. To verify the mass equation (|l.laj) . we decompose Rt and M x into two parts: an absolutely 
continuous measure and a singular measure. We distinguish between two cases. 

(1) . If y(x—, t) = y(x, t), one can trace x backward along the characteristic at least for a short 
time, and then applies the argument of Step #l(a) in the proof of the attractive case, k < 
case. 

(2) If y(x—,t) < y(x,t), then there is a 5-shock. The dynamics of the 5-shock for k > is a 
little different from the one for k < in that the shock may split or disappear. But before it 
splits or disappears, the argument of Step #l(b) in the proof of k < case is still valid. If the 
5-shock splits, we apply the argument to each sub-shock, the conclusion still holds. 

In a similar way, we can verify the moment equation for k > 0. □ 

Remark 3.1. Using the same technique, we can easily extend Theorem \l.l\ to measure density 
initial data, that is, Theorem \1.2l 

Remark 3.2. There is a global smooth solution if and only if x(a, t) = a + uo(a)t + ^KEo(a)t 2 
remains a monotonically increasing function with respect to a, that is, the solution remains 
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smooth as long as 



(3.2) = l + u' (a)t + -K Po (a)t 2 > 0, Vael,i>0. 

(j(y. 2 



Solving \3.2\l . we obtain the following critical condition: the repulsive Euler-Poisson system 
admits a global smooth solution if and only ifu' > — y/2npQ, [EnLTOi] . If the initial configuration 
fails to satisfy this critical threshold condition, then there are two critical times associated with 
the repulsive Euler-Poisson system. At the first critical time 

t Cl = inf{t|3a G R, such that 1 + u' (a)t + -Kpo(a)t 2 = 0}, 

the solution loses its C 1 smoothness; but after the second critical time, 

t C2 = inf{t|t > t Cl , 1 + u' (a)t + -Kp (a)t 2 > 0, Ma G R}, 

the solution given by hi. 6}) becomes C 1 again. 

Remark 3.3. The 5-shock splits for k > because the particles have "memory" back to time 
t = 0. More precisely, for k > 0, we have the following two constructions which may yield 
different solutions at time t: 

(1) Use hi. 6]) with initial data (po,uo) to compute p(-,t) and pu(-,t) directly. 

(2) Fix < t\ < t. First use \ 1.6]) with initial data (po,uq) to construct p(-,ti) and pu(-,ti). 
Then, solve the Euler-Poisson system subject to initial data p(-,ti), pu(-,t\) to obtain p(-,t) and 
pu(-,t). That is, apply 111.6]) with initial data (p,u)(-,ti) to obtain p(-,t) and pu(-,t). In this 
second construction, particles are losing their "history" before t±. That is, if particles collided 
before t\, they will "forget" they used to be separate and therefore stick together forever. 

To enforce the particles to lose their "history" at every moment (hence once collision happens, 
the particles will forget they were separate, and stick as one particle thereafter), one can impose 
the one sided Lipschiz condition of lemma \2.3l as an entropy condition: to this end, 111.5]) should 
be changed into 

(3.3) y(x,t) := sup ly y= arginf f Q x ,t{s)po{s)ds 

V I y£«\D(t) JO 

where D(t) = {z\3x e R,3t < t, s.t. y(x-,t) < z < y(x,t)}. 

4. Weak solution of multi-dimensional system with symmetry 

In this section, we extend the result to the multi-dimensional system with symmetry with 
explicit formulation of global weak solutions for the weighted Euler-Poisson system outlined in 
theorem 11.31 

Proof, of theorem 11.31 Since the initial data of (jl.9p are spherically symmetric, the solution of 
(jl.9p will remain spherically symmetric, p(x, t) = p(r,t), u(x, t) = w(r,t)*. Plugging these into 
(jl.9ap . we obtain 

u( ow] xx 1 x n 

Pt + J ^ r L {-,-) + pw{ — z-,x)+pw- = 0, 

that is 

, , . d(pw) n — 1 

(4.1) pt + + pw = 0. 

or r 

Multiplying (j4. 1 j) by r n_1 , we have 

Or 



(4.2) (r -v )l+ 3rrw = „. 
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8V X 

Since p is spherically symmetric, so does V, i.e., V(x, t) = V(r,t). Therefore, W = — , and 

or r 

AV = Vrr + Vr = p. Multiplying V rr + V r = p by r"" 1 , we obtain (r n -V r ) r = r n ~ l p. 

Plugging p(r,t), u(r,t) and V(r,t) into ()1.9b|) . we have 



pw-) + V • (pw 2 — x®x) = npr n ~ l VV, 



-) +V-(pw 2 \ 

r/t \ r z 

which is 

-{pw) t + ^-{n + l)x + (^-) -(x,x)x = Kpr^Vr-. 
r r z \ r z /rr r 

Further simplification yields 

(pw) t + P —(n + 1) + (^) r 2 = (pw) t + ^(n - 1) + (pw 2 ) r = npr n - l V r . 
r \ r z / r r 

Multiplying the above equation by r n_1 , we obtain 
(4.3) {r n ~ l pw)t + (r n -W)r = *r n ~ l pr n ~ l V r . 

Therefore, let f(|x|,i) = |x| n_1 / o(|x|, t), then (|4.ip and (|4.3p can be rewritten as 
q t + {sw) r = 0, 

(qw)t + (<V) r = K?r n V r =: K ?J E, — = i = ? . 

Consider (I4.4D with symmetric initial density r, 0) = <r(r, 0) = r n ~ 1 /o(r, 0) and anti-symmetric 
initial velocity w(—r,0) = —w(r,0). Let 

(4.5) E(r,t) = ^J^ s(a,t)da + J q(s,t)ds 

The difference between the electric field fj4.5|) and the previous one (which is defined by E(s, t) = 
\{§q s(s,t)ds + J Q r+ ?(s, t)ds)) is a constant f_ oo <^(s,0)ds. Physically, it corresponds to a 
Galilean transformation. For symmetric initial data, choosing E as (|4.5D is natural and conve- 
nience, since in such setting the particle located at the origin will not move and (?, w) will stay 
symmetric. Applying theorem 11.21 with the electric field (14. 5j) . we obtain a weak solution for 

(s,w). Then we recover ( p(x, £), (u(x, i) ) from mr, t), to(r, i) ) , r > 0. □ 



5. Concluding remarks 

We have constructed a global weak solution for the ID pressureless Euler-Poisson system. 
For the weighted multi-dimensional pressureless Euler-Poisson system with symmetry, which is 
essentially a ID system, we have constructed a global weak solution in the same manner. The 
open question is: is it possible to at least extend the method to the real multi-dimensional Euler 
system and obtain a weak solution from the physical laws directly? 
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